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Abstract. We consider the bi-Hamiltonian representation of the two-component coup- 
led KdV equations discovered by Drinfcl'd and Sokolov and rediscovered by Sakovich and 
Foursov. Connection of this equation with the supersymmetric Kadomtsev-Pctviashvilli- 
Radul--Manin hierarchy is presented. For this new supersymmetric equation the Lax repre- 
sentation and odd Hamiltonian structure is given. 
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1 Introduction 

The scalar KdV equation admits various generalizations to the multifield case and have been 
often considered in the literature [H [21 [3l HI [5]. However, the present classification of such 
systems is not complete and depends on the assumption which we made on the very beginning. 
Svinolupov [U has introduced the class of equations 

ui = uLx + a'i,k'^^'^x, (1) 

where i,j,k = 1,2, . . . , N , m are functions depending on the variables x and t and a* ^ are 
constants. The a*- ^ satisfy the same relations as the structural constants of Jordan algebra 

ai,fc(a^,ram,s - am.rOn.s) + cyclic(j, k, m) = 0. 

These equations possess infinitely many higher generalized symmetries. 

Giirses and Karasu [2] extended the Svinolupov construction to the system of equations 

Uf = bjU^^^ -\- Sj i^u-^u^, (2) 

where 6*, s*- ^ are constants. 

In general, existence of infinitely many conserved quantities is admitted as the definition 
of integrability. This implies existence of infinitely many generalized symmetries. Giirses and 
Karasu, in order to check the integrability of the system of equations ([5]), assumed that the 
system is integrable if it admits a recursion operator. Assuming the general form of the second 
and fourth order recursion operator they found the conditions on the coefficients bp s* ^ so that 
the equations in ([2]) are integrable. 

*This paper is a contribution to the Proceedings of the XVIIIth International Colloquium on Integrable Sys- 
tems and Quantum Symmetries (June 18-20, 2009, Prague, Czech Republic). The full collection is available at 
|http://www.emis.de/journals/SIGMA/ISQS2009.html| 
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A quite different generalization of multicomponent KdV system has been found by Antono- 
wicz and Fordy [3] considering the energy dependent Schrodinger operator. Ma [1] also presented 
a multicomponent KdV system considering decomposable hereditary operators. 

Several years ago Foursov [5] found the conditions on the coefficients 6*, s* ^ under which the 
two-component system ([2]) possesses at least 5 generalized symmetries and conserved quantities. 
He carried out computer algebra computations and found that there are five such systems which 
are not symmetrical and not triangular. Three of them are known to be integrable, while two of 
them are new. Foursov conjectured that these two new systems should be integrable. However, 
it appeared that one of these new systems is not new and has been known for many years. 
Drinfel'd and Sokolov in 1981 [6] presented the Lax pair for one of these new equation and hence 
this equation is integrable. 

In this paper we present the bi-Hamiltonian formulation and recursion operator for the new 
equation. These results have been obtained during the study of the so called supersymmetric 
Manin-Radul hierarchy. The application of the supersymmetry to the construction of new 
integrable systems appeared almost in parallel to the use of this symmetry in the quantum field 
theory. The quantum field theories with exact correspondences between bosonic and fermionic 
helicity states are not the only basic ingredients for superstring theories, but have been utilised 
both in theoretical and experimental research in particle physics. The first results, concerned the 
construction of classical field theories with fermionic and bosonic fields depending on time and 
one space variable, can be found in [71 |8l [9l [10]. In many cases, addition of fermion fields does 
not guarantee that the final theory becomes supersymmetric invariant. Therefore this method 
was named as a fermionic extension in order to distinguish it from the fully supersymmetrical 
method which was developed later [11^ \T2\ [T3| [14} [TS] . There are many recipes how some classical 
models could be embedded in fully supersymmetric superspace. The main idea is simple: in 
order to get such generalization we should construct a supermultiplet containing the classical 
functions. It means that we have to add to a system of k bosonic equations kN fermions and 
A;(A^ — 1) bosons {k = 1,2,..., N = 1,2,...) in such a way that they create superfields. Now 
working with this supermultiplet we can step by step apply integrable Hamiltonians methods to 
our considerations depending on what we would like to construct. 

Manin and Radul in 1985 [lOj . introduced a new system of equations for an infinite set of 
even and odd functions, depending on an even-odd pair of space variables and even-odd times. 
This system of equations now called the Manin-Radul supersymmetric Kadomtsev-Petviashvili 
hierarchy (MR-SKP). It appeared that this hierarchy contains the supersymmetric generalization 
of the Korteweg-de Vries equation, the Sawada-Kotera equation and as we show in this paper 
two-component coupled KdV equations discovered by Drinfel'd-Sokolov. 

2 Two- component KdV systems 

Let us consider a system of two equations 

ut = F[u,v], vt = G[u,v], 

where F[u, v] = F{u, v, Ux,Vx, ■ ■ ■) denotes a differential polynomial function of u and v. 

By the triangular system we understand such system which involves either an equation de- 
pending only on u or an equation depending only on v while by the symmetrical we understand 
such system in which G[n, f ] = F[v,u]. 

Definition 1. A system of t-independent evolution equations 



ut = Qi[u,v], vt = Q2[u,v] 
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is said to be a generalized symmetry of ([T]) if their flows formally commute 



Dk(Q)-Dq(K) = 0. 



Here Q = iQi,Q2), K[n, = {F[u, v], G[u, v]) , and Dk denotes the Frechet derivative. 

The first three systems in the Foursov classification are known to be integrable equations 



and are 




Ut 


= Uxxx + QuUx — 12vVx, 


Vt 


— '^^XXX ^UjVx^ 


Ut 


~ Uxxx ~l~ 3uUx ~\~ 3vVx, 


Vt 


= UxV + UVx] 


Ut 


= Uxxx + 2un^ + UVx, 


Vt 


= UUx- 



The first pair of equations is the Hirota-Satsuma system [16], second is the Ito system |17] . 
third is the rescaled Drinfel'd-Sokolov equation [5]. 

The fourth system of equations is a new one founded by Foursov 

Ut — Uxxx ~t~ Uxxx ~t~ '2iVUx ~t~ '^UVx, 

Ut = Uxxx — 9uux + 6vux + 3uvx + 2vvx- 

Foursov showed that this system possesses generalized symmetries of weights 7, 9, 11, 13, 15, 17 
and 19, as well as conserved densities of weights 2, 4, 6, 8, 10, 12 and 14, and conjectured that 
this system is integrable and should possess infinitely many generalized symmetries. 
The last system in this classification is 

Ut = 'iUxxx + 3fxxx + 4mMx + UUx + '^UVx, 

Vt = 3Uxxx + Uxxx — ^VUx — 2uVx — 2vVx, (3) 



and has been first considered many years ago by Drinfel'd and Sokolov [6] and rediscovered by 
S.Yu. Sakovich [l8]. 

Let us notice that the integrable Hirota-Satsuma equation has the following Lax representa- 
tion [19] 



while the integrable Drinfel'd-Sokolov equation possesses the following Lax representation [21 [20] 



L={d^ + {u- v)d + (n, - Vx)/2) + {u + v)d + {ux + Vx)/2) , — = 4 [1%^ L] . 



On the other side, the Lax operators of the Hirota-Satsuma equation and of the Drinfel'd- 
Sokolov equation could be considered as special reduced Lax operators of the fourth and sixth 
order respectively. Indeed, the Hirota-Satsuma Lax operator could be rewritten as 



L = [d"^ + u + v){d'^ + U-V 



dL 

dt 



4[Lr,L] 



L = d^ + g2d^ + gid + go 



(4) 



where 



92 = 2u 



gi = 2{ux - Vx) 



go = Uxx + U^ - Uxx 



— v 



,2 
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In this context, one can ask what kind of the equations follows from the fifth-order Lax 
operator which is parametrised by two functions of same weight. Let us therefore consider the 
following Lax operator 

L = d^ + h2d^ + hsd"^ + h^d + h5, 

where hi, i = 2,3,4,5 are polynomials in u and v and their derivatives of the dimension i. 
Computing the Lax representation for this operator 

we obtained then 

L={d' + lud + in.) {d' - in) (5) 

produces the system of equation ^ . 

Let us notice that the Lax operator ([4]) is factorized as the product of two Lax operators. 
The first one is the Lax operator of the Kaup~Kupershmidt equation while the second is the 
Lax operator of the Korteweg-de Vries equation. It is exactly the same Lax operator which has 
been found by Drinfel'd and Sokolov [6]. 

Hence we encounter the situation in which the Lax operator of the Korteweg-de Vries and 
the Kaup-Kupershmidt equations can be used for construction of additional equations. This 
could be schematically presented as: 







Lkk 




Hirota-Satsuma 


equations ([3]) 


Lkk 


equations ^ 


Drinfel'd-Sokolov 



where LKdV, LKdv are two different Lax operators of the Korteweg-de Vries equation while Lkk 
and Lkk are two different Lax operators of the Kaup-Kupershmidt equation. 



3 The recursion operator and bi-Hamiltonian structure 

From the knowledge of the Lax operator for evolution equations one can infer a lot of properties 
of these equations. The generalized symmetries are obtained by computing the higher flow of 
the Lax representation while the conserved charges follow from the trace formula [21] of the Lax 
operator. 

Using this technique we found first three conserved quantities for the equation ([3]) 
Hi = j dx{v'^ +4u'^ +6uv), 

H2 = j dx(495u4xu - 510ulu + 32m^ + 2v'^ + 630v4xU + 180^4^.^ - 210v^^vu - 210v^.n 

= / ..(I82250„,.u + ™%00„.....„ + + 259200„L"=+ 223425„„„?. 

— 104400n^u^ + 1344n^ + 222750^8x1* + IQSlbvsxV - 148500n6xTO - im^lbv^^u^u 

— 594000^52^^ f.u — l\l?>7hQv^xVxxU — ^2S2hQv^xVxxV + 54450u5a;f — S2hv^xVV?' 

— 742500^^2,2, u- 74250f^2.2.u- QlSlbvxxxUxxxU— 2VlS0QvxxxUxu:^ + 2Q73(^{)VxxxVxVU 
+ 70125?;^ n2 _^ 19575^2 2 _^ i63350t;2 2^ - l^2,QbQvxxulu + 2^7mQvxxvlu 
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-'xx ^x 



32175< 



19800t;?u^n + 13200v? 



vu 



185625?;2;U5a;U 



+ ISSlOOvxUxxUxU — 79200vxUxU^ + 825vxv'^UxU + 57750vxVUx 
+ 198v^u + eeOu'^M^ - 7425?;^ u^^u + UOv^u^ + 44550t;2u4xii - UbbOv'^ulu 



llU7hvu(ixU - mfiQvulxU - 59400vu>^ + 3168um^). 



Taking into the account a simple form of the first Hamiltonian it is possible to guess the first 
Hamiltonian structure 



di 



P 



5u 
6Hi 



5v 



I 



j^Qu + ud 

3(9^ - 2(dv + vd) 




In order to define the second Hamiltonian structure we first found the recursion operator. We 
used the technique described in [22] and we found the following tenth-order recursion operator 



R 



-^^5^° -F268 terms 
-^51° 268 terms 



-^(9^°-F268 terms 
-^d^^ ^2<o8 terms 



Next we assumed that this operator could be factorized as R = J~^P where is the inverse 
Hamiltonian operator. Due to the diagonal form of the first Hamiltonian structure it is easy to 
carry out such procedure and as a result we obtained the second Hamiltonian structure 



dt 



SH4 

Su 
5v 



where 



i74 = J dx{21uiQxU + 26fioa;?/ -|- 95 terms) 



and the explicit form of H4 and J ^ is given in the appendix. 



4 The derivation of the Lax representation 

The Lax operator of equations ([5]) has been discovered accidentally during the investigations of 
the super symmetric Manin-Radul hierarchy. This hierarchy can be described by the supersym- 
metric Lax operator 

00 00 

L = V + fo + Y,bjd~^^ + Y, fjd~' , (6) 
j=i j=i 

where the coefficients bj, fj are bosonic and fermionic superfield functions, respectively. We shall 

use the following notation throughout the paper: d and T> = -\- 6d. As usual, (x, 9) denotes 

= 1 superspace coordinates. For any super pseudodifferential operator A = ^ ^jj'iP'' 

i 

subscripts it denote its purely differential part = ^ ^jj'iP^ or its purely pseudo-differential 

i>o 

part A^ = ^ a^ji2T>~^ respectively. For any A the super-residuum is defined as Res A = a-i/2- 
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The constrained (r, m) supersymmetric Manin-Radul hierarchy |25j is defined by the fol- 
lowing Lax operator 

1 m 
j=0 j=0 

This hierarchy for even r has been widely studied in the literature in contrast to the odd r which 
is less known. Further we will consider this hierarchy for odd r = 3, 5 and m = 0, T = \I/ = 0. 

The Lax operator for r = 3 and m = 0, T = ^' = has been considered recently by Tian 
and Liu 



L = p3 _^ $, 

where $ is a superfermion function. Let us consider the following tower of equations 

Lt,k = 9[L,L':/^]. 
The first four consistent nontrivial equations are 



' X' 



^t,10 = ^5x + ^^xxx^l + 5^xx^l,x + 5^x^l, 
^r,ll = $l,5x + 3^1,^^^$! + Q<^l,xx'^l,x + '2<^l,x'^i 

- 3$4x^> - 2<^xxx^x - 6^'xx$$l - 6^>:^^>^>i,^, 

where t is a usual time while r is an odd time. 

The third equation in the hierarchy in the component = ^ + 9w reads 

6 = ibx + 5wC XXX 

Wt = W5x + ^WWxxx + ^WxWxx + 5w'^Wx - ^^.xxx^x 

and it is a supersymmetric generalization of the Sawada-Kotera equation. This equation is a bi- 
Hamiltonian system with odd supersymmetric Poisson brackets [23]. The proper Hamiltonian 
operator which satisfies the Jacobi identity and generates the supersymmetric = 1 Sawada- 
Kotera equation is 

where Hi = f ^^x dxdO and 

P = (Vd^ + 2(9$ + 2$5 + V^V)d-'^ {Vd'^ + 2d^ + 2^d + V^V) . 
The implectic operator for this equation was defined in [23] as 

J$t = ^, J = dxx + {V^)-d-\V^)x + d-^<^xV + ^xd-^V, (7) 

where 
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This supersymmetric equation possesses an infinite number of conserved charges [23] which are 
generated by the supertrace formula of the Lax operator. However, these charges are not reduced 
to the known conserved charges in the bosonic hmit. Hence we can not in general conclude that 
from the supersymmetric integrability follow the integrability of the bosonic sector. 
Let us now consider the Lax operator ([6]) for r = 5 and m = 



where U and V are superfermionic functions U = ^ + 6u, V = i{j + 9v. The first nontrivial 
equations in the hierarchy generated by this Lax operator is given as 



Ut = Wxxx + ^Vxxx - 2Ux{VU + VV) + U{6VUx + - V^W + V{3VUx + W,), 

Vt = 3Uxxx + Vxxx + SUxVU - U{SVUx - 6VVx) + Vx{iVU + VV) - V{iVUx + 3PK). 



The bosonic sector of the latter system where ^ = 0, "0 = gives us the system of two interacted 
KdV type equations discovered by Drinfel'd-Sokolov. 

Interestingly, the Lax operator equations ([7]) did not reduce in the bosonic sector to our Lax 
operator ([5]), however, its second power reduces that one can easy verify. As we checked, this 
system possesses the same properties as the supersymmetric Sawada-Kotera equation. Namely, 
this model, due to the Lax representation, has an infinite number of conserved quantities, which 
are not reduced to the usual conserved charges in the bosonic limit. For example, the first two 
conserved charges are 



We found the following odd Hamiltonian structure for our supersymmetric equation ([8]) 



Unfortunately, we have been not able to found second Hamiltonian structure for our superequa- 
tion. 

A Appendix 

The conserved quantity i/4 is 



L = p5 + l{dU + Ud) - \VVV 




J dxde{75UxxxU + 32UxU{VU) - 24f/^C/(Py) + dOV^xxU + SOV^xxV 
+ 36VxU{VU) - 6VxU{VV) - AVxV{VV) - ?>WU{VVx)). 



(8) 





~l~ 315 VxxxU5xU~\~ '^ggVxxxVxxUxU-\- g3 VxxxVxx'^x'U ''^xxx^xx^x^ ^j2^^^^^^^^ ^ 
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^qy'^'xxx'^'x'^^U g^r^VxxxVUxxxU -\- 27 "^xx"^ gi^'^xx'^ g45'^xx'^xx'U' ~^ i^if^V^^U 

- _88_,,2 3 _ J592. 2 2 _ 26_ 2 2 , 2782 „ ^, , 2548,, n - '^V v'^VU 

2025 xx^ 14175 ^xx*^ " 135 '^xx'^" ' 45 ^xx^exU-T^ 135 ^JxxU-xx"- 2025^^^^x'^ 

_ 2236,, ,, ,, ,, I 10£,, ,,4,, I ^08_ 3 2 , 104 2 _ J9£,;4,, 

2g35 i^xx L'x ""xxx 6075 ^^ 42525 a:x'^ " T ^Qgt'xx" "-xx"- 2025 x" 

^ 6075 a; 567 ^x "x 4725 ^x^xx" 1575 ^x'^x" 14175'^x"' 6075 a; 

+ lis^x^^^^ - ^^x^^^^^ + if ^^x^^^xxW - ^m'^l'"^^ + ^^^x'W7xU - ^VxUAxUxU 

945 ^x^xxx^xxU "4725" ^x'U 2^1* ~l~ 4725 '^xVU^xU 14175 VxVUxxUxU 3189375 ^ 
104 „,6„, 416 „,5„,2 I 104 „,4„, „, 832 „,4„,3 1144 „,3„, „, , 1664 „,3„, „ ,2 



455625 



^ 455625^ ^ + 8505^ '"'^^^ 637875^ ^ 14175"" W.4xtt + 42525^ '"'ajx 

, ^16_ 3 2 , 1664 3 4 , 338 2 , ^7,2,, , ,,2 1 101-^2,, ,, ,, 

14175 x" ^ 637875 945 ""Dx^t^ 4725 "'4x"' 1 4725 "'xxx^x"' 

_ i01,;2,,2 , 1661 2 3 , 10816 2 2 2 , 3328 2 5 _ 247,,,,„ ,, + 208 : 
4725 ""xx^-T^ 14175" "'xx" 42525 a;" ^ 455625 33^5 f^"'8x " ]^o5 6x" 

_ 10868 „, _ 14872 , 4576„,„. „,3 _ 29692 2 , 102128 „.,, ,,2 

945 ""'5x"'X"' 315 "■4x "'XX 4725 ""'4x"' 945 ""'xxx"' 14175 '^"'xxx"'x"' 

+ 65416,,,,2 2 . 123682 ,,2 , 1664,,,, .A 1 4576 2 3 , 3328 6 
^ 14175 '^"■xx"' 14175 'J'J'xxU'x"' ' 5975 '^"-xx"- 6075 a;" 455625" ■ 



The inverse Hamiltonian operator J ^ has the following form 



( 

-I _ I "^1,1 "^1,2 

\ "("^1,2) ^^2,2 



where 

•^iTi = - lil^^ + «i,i,59^ + ai,i,35^ + ai,i,i9 + 61,19-^ + 61,143-^61,1,2 - h.c, 
5 

Jl,2 = -ifE^' + Yl "l'2,i5' + 61,25-1 + d-^Ci,2 + 6l,2,l5-l6i,2,2, 
i=0 

-^2,2 = -tEo^'' + 0-2,2,58^ + a2,2,3d^ + a2,2,id + 62,29-1 + 62,2,15-^62,2,2 - h.c, 
01,1,5 = (-43'u + 14i;)/1125, 

01.1.3 = (2251*^^ - 448^2 - 690vxx - TTd^ - 42to)/16875, 

01.1.1 = ( - 3150n4x - ISOOuxxU + 2397n^ - 888n^ + 720?;4x - 2520vxxU + ^20vxxV 

+ 7987;^ - 1512vxUx + 56v^ - 126v\ - 728t;n2) /101250, 
61 1 = ( - 6750u6x - 9540u4xU - 19080tta:xx'«x - 14310'uL - 5520uxxU^ - 5520n^ti - 352n^ 

- 4050vgx - 3780t;4x^t + 1890?;4xt; - 8107;3j:Ua; + 6480i;3a;Ua; + 3510t;^^ + 2160?;a;xiixx 

- 21Q0vxxu'^ - 540vxxv'^- liAOvxxVU+ 720vlu - 810't;^i;+ 59i0vxUxxx - WSOvxUxU 

- 2160vxVUx + 18t;^ + 96v^u - WSOv^Uxx - 288v^u^ + 2970^4^; - WSOvUxxU 

- 540vul - 576vu^) /759375, 

6i,i,i5"i6i,i,2 = ( - 990'U4x - I020uxxu - 510ul - 128n^ - 630t;4x - 420vxxU 
+ 210vxxV + 210vl + 210vxUx - Uv^ - 56v^u + 210t;Ua;x 
- 168^2) (4n + 3w)/759375, 

01,2,5 = (-38^ + 19t;)/1125, 

01.2.4 = (-36ti + 68u)^/1125, 

01,2,3 = {—4:5uxx — 124n^ + 515^2:2: — 31t;^ — 36vu) /5625, 

01.2.2 = (885u32: — 656uxU + 1680t;3a; + 178?;a;ti — 354^^;^ — 332to2:)/16875, 
01,2,1 = (2520'U4x - 1428ita:xii - 15991*^ - 328u^ + 2610d4x - l56vxxU 

- 957vxxV - 876vl - 1716f^Ua; + 41v^ + lUv'^u - mivuxx - 228vu'^) /50625, 
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fli,2,o = (l575u5a; + 534u3xU + '2'76uxxUx — 128uxu'^ + lOSOv^x + WSv^xU 

— 534v3xV — 873VxxUx — ^'^'OQVxxVx — i^SSVxUxx + SAVxU^ 

+ 96vxV^ + 166vxvu + I28v^ux - 792t;n3x + Svuxu) /50625, 
bi,2 = (3375^60; + 4770u4a;n + dSAOu^xUx + 7155u^^ + 2760«a;a;ti^ 

+ 2760ulu + 176^'' + 2025^6^ + 1890v4xn - 9457;4x't^ + 405u3a:'Ux - 3240v3xVx 

- 1755t)^2. — lOSOvxxUxx + 10807;a;a;ti^ + 270t'2,u'^^^ + 720vxxVU — 360v'^u 
+ A05vlv - 2970ua;ti3x + 540t;xUa;U + lOSOvxVUx - 9u^ - 4.8v^u 

+ 540u2ua;a; + 144?;2^2 _ 14851,1^4^ + 540vUxxU + 270vul + 288^^2) /759375, 
c^2 = (-4050n6x - 3780n4xii - USlOu^xUx - 90451*^ - 2160uxxU^ - 3780n^n - 144n^ 

- 2700t;6x- 2160?;4x'« + 1530v4xV- 4320v3xUx+ 3060^3^?;^ + 2295^^3,- lOSOvxxUxx 

- 2\mvxxU^ - 'i'^^vxxv^ - '^^'dvxxvu - 270vlu - 390vlv + WSOvxUsx - 4320vxUxU 

- 5A0vxVUx + llv'^ + 72v^u - bAOv'^Uxx + UAv'^u^ + 1890t;n4j: - IMOvUxxU 

- 360vul - 192vu^) /759375, 

fei,2,i9~^fei,2,2 = ((-990^4^ - 1020uxxU - bWul - 128u^ - 630^4^ - 420vxxU + 210?;^^?; 

+ 210';;^ + 2WvxUx + Uv^ - 56v^u + 210vuxx - 168vu^) {3u + v) 

- (4n + 3v)d"^ (630n4x + 420'u^^n + 525n^ + 56n^ + 360t;4x - '^50vxxV 

- 75vl + 8v^ + 42v'^u - 2WvUxx + 56vu'^)) /759375, 
02,2,5 = (-9n + 7w)/1125, 

02,2,3 = (540na;a; - 122^2 - 30vxx - 33z;^ - 28t;n) /16875, 

02,2,1 = (-1170'u4x + 9001*^^^ + 1203^^ - 120^^ + 5A0v4x - l^OvxxV 

+ \<o2vl - \'ii8vxUx + 30w^ + llOv^ti - 'ol<iivUxx + 120^^) /101250, 
62,2 = (4050u6x + 3780u4^u + 14310n3^ti^ + 9045u^^ + 2\mUxxU^ + 3780'u^w 
+ 144w^ + 27mv^x + 2\<o<^V4xU - V^Z'iiv^xV + 4320^3^ u^: - 3060^3^ -y^^ 

- 2295^;^^ + X^mvxxUxx + 2160^^^^^ + Z^VxxV^ + hmvxxvu + 270v^.n + 390u^z; 

- 1080fa;M3a; + 4320vxUxU + 540^2,^2: - llv'^ - 72v^u + M^v^Uxx 

- lUv'^u^ - 1890t;n4^ + lUOvUxxU + 360vul + 1927;^^) /1518750, 

ft2,2,lC?~^&2,2,2 = ( — 630tt4a; — 'i20UxxU — 525li^ — 56u^ — 360v4x 

+ I50vxxv + 75vl - 8v^ - \2v^u + 2\^vuxx - h<ovv?)d~^ {Zu + z;)/759375. 
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